Stable droplets and growth laws close to the modulational instability of a domain wall 
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We consider the curvature driven dynamics of a domain wall separating two equivalent states in 
systems displaying a modulational instability of a flat front. We derive an amplitude equation for the 
dynamics of the curvature close to the bifurcation point from growing to schrinking circular droplets. 
We predict the existence of stable droplets with a radius R that diverges at the bifurcation point, where 
a curvature driven growth law R(t) ~ t 1 ^ 4 is obtained. Our general analytical predictions, which are 
valid for a wide variety of systems including models of nonlinear optical cavities and reaction-diffusion 
systems, are illustrated in the parametrically driven complex Ginzburg-Landau equation. 
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Growth of spatial domains of different phases in the 
transient regime of a system approaching thermody- 
namic equilibrium has been studied long ago [|lj|. This 
phenomenon is a prototype of nonlinear dynamics of 
an extended system governed by the motion of domain 
walls. Power laws for the growth of a characteristic size 
and self-similar evolution (dynamical scaling) have been 
established. Clear physical mechanisms explaining the 
different asymptotic growth laws have also been identi- 
fied [|l]]. It is convenient to consider domain growth in 
systems with conserved and non-conserved order param- 
eters separately. In the first case one talks of spinodal de- 
composition leading to the Lifshits-Slyozov t 1 / 3 power 
law. In the case of non-conserved order parameters, the 
dominant mechanism is curvature driven minimization 
of surface tension, leading to the Allen-Cahn (AC) t 1 / 2 
power law [^] . 

The study of domain growth and domain wall motion 
in non-equilibrium systems that do not approach thermo- 
dynamic equilibrium is by far more complicated. There 
is presently only a partial understanding of a variety of 
possible situations [|3jQ]. In particular a number of re- 
sults have been recently reported for domain growth and 
domain wall motion in the transverse plane of nonlinear 
optical systems However, these numerical results 

often fail to identify the physical mechanism responsible 
for the observed dynamics. Moreover, the power laws 
for domain growth do not always correspond unambigu- 
ously to an asymptotic and self-similar dynamics. An 
exception is the i 1 / 2 law obtained for two different opti- 
cal systems for which evidence of dynamical scaling has 
been given . Visual inspection of the time dependent 
configurations clearly indicate curvature driven domain 
wall dynamics, but no derivation of such a law has been 
given for these systems. Surface tension, whose mini- 
mization may be used to explain the AC law, is not an 
appropriate concept for these optical systems. 

In this letter we analyze the transition from a regime 



characterized by a t 1 ' 2 growth law to one of labyrinthine 
pattern formation due to a modulational instability of 
a flat domain wall connecting two equivalent homoge- 
neous solutions. Such a transition from coarsening to 
labyrinthine regimes has been observed experimentally 
in reaction-diffusion jl0| ] and in optical 1 11 ] systems, and 
numerically in [g], as well as in Swift-Hohenberg models 
Jl2|]. At difference with systems approaching thermody- 
namic equilibrium, the coefficient of the t 1 / 2 power law 
can now change sign upon variations of a control param- 
eter p at p — p c - For p < p c ( p > p c ) a circular domain 
of radius R of one solution embedded in the other grows 
(shrinks as in the AC regime) as R(t) = y/R(Q) — -ft 
with 7 < (7 > 0). Moreover, for p < p c a flat front 
is modulationally unstable and a generic initial condition 
develops in a labyrinthine pattern. Close to the modu- 
lational instability at p c we derive an amplitude equation 
for the curvature of circular domains. At p c a circular do- 
main grows as R(t) ~ t 1 / 4 . For p larger but close to p c 
an initially small (very large) domain grows (shrinks) un- 
til a stable droplet (SD) is formed. The radius of the SD 
diverges at p c . The SD should not be confused with other 
types of soliton-like localized structures (LS) appearing 
in the AC regime when the front oscillatory tail prevents 
the domain collapse |Q,|l3j^]. The SD are also different 
from the critical nucleus of nucleation theory |Jl|] where 
the droplet is unstable and inner and outer solutions are 
not equivalent. Note also that droplets found in reaction- 
difussion systems [ |l4| | involve two non equivalent states. 

Our analytical results are illustrated by considering a 
prototypical model, the parametrically driven complex 
Ginzburg-Landau equation (PCGLE) [ |L5| , |l6[ l, which de- 
scribes a oscillatory system parametrically forced at 
twice its natural frequency [|l0|] : 

d t A = (1 + ia)V 2 A+ + iv)A (1) 
-{l + i(i)\A\ 2 A + pA* , 

where p measures the distance from the oscillatory insta- 



1 



bility threshold, v is the detuning parameter and p > 
is the forcing amplitude. For p ~ v ~ a large compared 
to other parameters, a finite wavelength instability form- 
ing hexagons takes place for p < ph, while for p > ph 
there are two equivalent stable homogeneous solutions 
(frequency locked solutions) [jig]. 

We consider a real N component vector field \5(x) 
whose dynamical evolution in two spatial dimensions is 



d t V = DV 2 * + W(*,p) 



(2) 



where the matrix D describes the spatial coupling, W is 
a local nonlinear function of the fields and p a control 
parameter. Eq. (Q) is invariant under translations and un- 
der the change x — ► —x. In addition we assume that 
it has a discrete symmetry Z that allows for the exis- 
tence of two, and only two, equivalent stable homoge- 
neous solutions, and that, in a Id system, they are con- 
nected by a stable Ising front ^o(x,p). An Ising front 
satisfies ^0(^0 — x) = Z^ (x — xq) where Xq is the 



front location [17], thus the Id front (and equivalently a 
flat front in 2d) is stationary, £>V 2 * + W($o,p) = 
0. For the PCGLE = (Re[A(x)}, Im[A(x)}), 

D = ((1, a) T , (-a, 1) T )), $o(x) describes a flat front 
connecting the two homogeneous solutions (we consider 
p > ph) and Z — —I where / is the identity matrix. 

Let X(s,t) represent the instantaneous position vec- 
tor of the front in the x plane, where s is the arclength. 
It is convenient to define a coordinate system (r, s) that 
moves with the front such that x = X(s, t) + rf(s, t), 
where f is a unit vector normal to the curve X and the 
coordinate r is the distance of the point x to the front 
In the moving frame Eq. (ph becomes 
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(3) 



where v = dtX ■ f is the (normal) front velocity and 
9 = ks is the azimuthal angle. We analyze the dynam- 
ics of slightly curved fronts as a perturbation of the flat 
front #(r, s,t) = $o(r) + $i(r,s,t). We assume that 
i) kw -C 1, with k — V • f the curvature and w the 
front width, ii) in the moving frame the front profile de- 
pends at most weakly on t (|9t^| <C \nDd r ^>\) and iii) 
k is a function which depends at most weakly on s, thus 
|«;<9g\&| -€i |9r$|. Linearizing around ^0 we have 



-(vl + KD)d r $ , 



(4) 



where M\ = D l M + 5^,jW% . Due to the transla- 
tional invariance of (Eh M is singular, Meo — where 
eo = d r ^o is the Goldstone mode. The solvability con- 
dition applied to (0) leads to v = — 7(p)«, where 



1 r 00 

7(P) = f / oo • -D 
^ J —00 



e cir , 



(5) 



T = ao ■ eodr [ p^| ] and Sq is the null mode of Aft. 
For a circular domain k = 1/R and 



= R = - 7 (p)/R. 



(6) 



From Eq. (j4j) one obtains that the perturbation of the 
front, t) = K(t)<pi(r), is independent of s while 

the dependence on t comes only through the curvature. 
01 satisfies 



M0x 



-(-jI+D)e . 



(7) 



For systems such that the diffusion matrix is propor- 
tional to the identity, D = dl, 7 takes the constant value 
d independently of the profile of the front and any sys- 
tem parameter. In this case, the rhs of (j7h vanishes and 
if 1 must be either zero or proportional to the Goldstone 
mode eo ■ Physically this means that the fronts translate 
without changing their radial profile. The front veloc- 
ity is proportional to the curvature with opposite sign 
v = —dn, which is the well known AC law [Jl|^]. This 
law implies a coarsening regime with a i 1 / 2 growth law 
and shrinking of circular domains. 

In general, D = dl + C with a non zero matrix C. 
C leads to a contribution to 7 that depends on the pro- 
file of the front and therefore, on the system param- 
eters. From Eq. (Q) we have that (pi is no longer pro- 
portional to the Goldstone mode. This means that the 
transverse profile of the front is now deformed. Since 
*&i(r,t) = K(t)<pi(r), the amount of deformation is pro- 
portional to the curvature and C. The d contribution 
is generally positive and, for wide parameter regions, 7 
is also positive. Here flat walls are stable and the AC 
law still applies [^,0]. Circular domains shrink but the 
presence of oscillatory tails in the front may prevent the 
droplet to disappear forming a LS as found in nonlinear 
optical cavities [|7[pj|,^|| . 

The crucial point is that for some parameter values the 
contribution to 7 due to C may be negative and larger 
than d. 7 changes sign and a bifurcation occurs. This is 
particularly relevant in nonlinear optical systems where 
the spatial coupling is diffractive and d — 0. For 7 < 
the velocity has the same sign of the curvature leading 
to the growth of any perturbation of the flat wall. The 
condition for vanishing 7, f_ ao • De^dr = 0, is in 
fact the criterion for the modulational instability of a 
flat front. Modulational instabilities in fronts connecting 
two equivalent homogeneous states have been observed 
in different systems Starting from a random ini- 

tial condition, the system develops labyrinthine patterns 
fl[Io|]. Also, a circular domain of radius R grows like (^|) 
until its boundary breaks up because of the modulational 
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instability leading to the formation of a labyrinthine pat- 
tern. 

In Figjj] we show the value of 7 for the PCGLE cal- 
culated using definition (||), with the profile of the Id 
front calculated numerically by solving the stationary 
form d t A = in Id of (|J). This calculation identifies 
p c for which 7 = 0. 

We are now ready to obtain an amplitude equation 
for the curvature in the vicinity of p c . We start con- 
sidering the case of a circular domain wall (for which 
v = —k/k 2 ). Close to p c (7 w 0) we perform a mul- 
tiple scale analysis in e of Eq. (O) with p = p c + epi, 
* = *o + e 1/2 *i + e*2 + e 3 ^*3, K = e 1/2 Ki and 
dt = e 2 dr- At order e 1 / 2 we obtain = — KiDe . 

The solvability condition implies ao ■ De^dr — 0, 
which is automatically satisfied at p = p c . One finds 
"fi = Ki0i where ipi satisfies M(pi = —Dcq. Due to 
the symmetry of the front the solvability condition at or- 
der e is always fulfilled and we get <3>2 = Pi<f2 + nf03 
withM(^ 2 = -d p W\o mdMjtpl = -I>J(d r <pj-rej)- 
\8y W 1 \o(p{(p'l where | means evaluated at ^ and 
p c . At order e 3 / 2 , the amplitude equations for the curva- 
ture Ki and the radius of the circular droplet 

d TKl 3 

5-= C1P1K1 + (8) 

ft 2 

d t R = -ci(p-p c )— - c 3 -^3 (9) 

where 

1 f°° 

c\= p / a 0l {D)d r Lp 3 2 + S^jdpW'loipl + 

^ J —OO 

S&vkW^ovi^dr (10) 

1 r°° 

c 3 = = / a w [D , j (a r tp' a - r8r(p{ + r 2 e J ) + 
SyjykW^cKplip^ + ^8 9j9 k 9 iW i \o(p{(p'lip l 1 ]dr (11) 

are obtained from the solvability condition. The coeffi- 
cient ci is positive since 7 = c\(p — p c ), and we are 
considering 7 > for p > p c . If C3 is negative (super- 
critical bifurcation) our analysis predicts just above the 
modulational instability the existence of stable station- 
ary circular domains (SD) with a very large radius 

Ro = * (12) 

VP - Pc V ci 

In Fig. ^| (left) we show the form of the SD for the 
PCGLE. At the center the field takes the value of one 
of the homogeneous solutions, so the wall of this struc- 
ture is a heteroclinic orbit between the two homogeneous 
states. The radius of the SD can be extremely large di- 
verging at p = p c . Fig. H (right) displays the radius of 



the SD and that of the LS calculated solving numerically 
(l+ia)(d 2 +±d r )A+{n+iis)A-(l+il3)\A\ 2 A+pA* = 
0. The inset in Fig. ^] (right) shows the linear dependence 
of 1/Rq with p as predicted by ([L?]). In spite of the fact 
that there is a smooth transition from LS to SD, they are 
intrinsically different. The interaction of the oscillatory 
tails which is responsible for the existence of LS does 
not play any role in the SD. The stabilization mechanism 
comes from the counterbalance between the i?~ 3 contri- 
bution to the front velocity and the shrinking due to the 
R^ 1 contribution. 

At the bifurcation point p = p c , Eq. ^ becomes 
d t R = -c 3 /R 3 and any circular domain of one solu- 
tion embedded in the other grows as R(t) ~ t 1 / 4 . In 
Fig. |]we show the time evolution of the radius of a cir- 
cular domain for the PCGLE at p = p c . It fits nicely 
the theoretical dependence we predict. Close to p c , in 
the regime of existence of the SD, there is no asymptotic 
power law of domain growth since at very long times the 
SD is formed stopping the growth process. During the 
transient, an initially small (very large) circular domain 
will grow (shrink) following (^). 

So far we have considered the dynamics of domains 
with radial symmetry. When the system evolves from 
random initial conditions other dynamical mechanisms 
come into play. Close to the bifurcation point the 
main non-radially symmetric contribution to the veloc- 
ity comes from the variation of the curvature along the 
front. The derivation of the amplitude equation for the 
curvature leads to an additional term in Eq. d8j) 
where C2 = jt do ■ D<p±dr. The front velocity be- 
comes v = —ci(p—p c )K—C2K 2 d 2 K—C3K 3 . Consistently 
with our approximations, dgK will change at most at or- 
der k°, so the non-radial contribution is at least of order 
k 2 . Close to the bifurcation point the term proportional 
to k can be neglected, and the term given by — C2« 2 9 2 k 
might be dominant compared with the one proportional 
to K 3 . In this case the front velocity is proportional to 
ft 2 . However, for any closed boundary 2 k is positive 
in some parts of the wall and negative in others, so this 
term does not lead to an asymptotic growth law. For 
C2 < 0, which is the case for the PCGLE, this term tends 
to reduce the curvature differences, so at p c an arbitrarily 
shaped domain first becomes circular until the contribu- 
tion of d 2 K vanishes and then the circular domain grows 
as R(t) ~ i 1 / 4 due to the C3 term (see Fig. Q). 

In summary, we have analyzed a generic situation of 
domain wall motion driven by curvature effects in which 
the proportionality coefficient 7 between wall velocity 
and curvature changes sign at a bifurcation point. In op- 
tical systems, this change is a consequence of the diffrac- 
tive coupling between real and imaginary parts of the 
complex field amplitude. The amplitude equation for 
the curvature in the vicinity of the bifurcation point pre- 
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diets the existence of stable nonlinear solutions which are 
droplets of one phase embedded in a background of the 
second equivalent phase. Nonlinear dynamics of the cur- 
vature leads to growth laws different from the AC in- 
growth law. The existence of a large characteristic length 
given by the radius of the SD destroys the possibility of 
self-similar evolution. 
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FIG. 3. Growth of a circular domain as function of the time 
atp — p c - The symbols correspond to the numerical integration 
of Eq. (|l|), the solid line is the theoretical prediction (g) with 
C3 = 0.3129 calculated from ( |ll| ) and the dotted line has a 
slope 1/4 showing the asymptotic behavior. 




FIG. 4. Evolution of an arbitrarily shaped domain atp = p c . 
From left to right: t=0, t=2400, t=80000 and t= 1650000. 
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FIG. 1 . Value of the 7 as a function of the forcing amplitude 
p for the PCGLE calculated from (j^). We take here a — 2, 
/3 = 0, v = 2 and n = 0. For these values of the parameters 
p h = 2.09 and p c = 2.56629. 
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FIG. 2. Left: Spatial dependence and transverse section of 
a SD for the PCGLE. Right: Radius of SD (solid line) and LS 
(dotted line) as a function of the forcing for the PCGLE. The 
inset shows the linear dependence of 1 /Rq with p close to the 
bifurcation point as predicted by Ql 2p . 
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